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Abstract 

We have calculated the energies of the phase-separated states for 
degenerate antiferromagnetic semiconductors including the possibility 
of the existence of conduction electrons in the antiferromagnetic part 
of the phase-separated states. It is demonstrated that, at T = 0, 
the minimum energy corresponds to a droplet phase with absence of 
electrons in the antiferromagnetic part. 

1 Introduction 

Degenerate antiferromagnetic semiconductors are obtained by strongly dop- 
ing antiferromagnetic semiconductors (e.g., EuTe or lanthanium mangan- 
ites). Over a concentration range of doping impurities, the ground state of 
these compounds will be a mixture of antiferromagnetic (AF) and ferromag- 
netic (FM) phases |], @, f|. Starting from a pure compound and doping, it is 
found that the ground state is AF and non-conducting up to some conduction 
electron concentration, driven by impurities, n^, for which the homogeneous 
state turns to be unstable against phase separation. The ground state be- 
comes then inhomogeneous with a simply connected AF phase and a multiply 
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connected FM phase. The exact geometry of the multiply connected phase 
is of fractal nature. On increasing the doping, a geometric transition takes 
place in which the topology of the sample changes due to percolation of the 
FM phase. Now, the ground state corresponds to a multiply connected AF 
phase and a simply connected FM one. The doping concentration at which 
this transition occurs is denoted by At this point an important change in 
the conductivity is expected, because the material becomes a conductor. On 
a further increase of doping, a certain concentration, n u , exists at which this 
phase-separated state starts to be unstable and the sample becomes again 
homogeneous but now a FM. Transitions between these phases are also ob- 
served by varying magnetic field and temperature [|J. 

The physical reason for having a phase-separated state in an antiferro- 
magnetic semiconductor is the dependence of the energy of the charge carriers 
that appear on doping on the magnetic order of the lattice. The charge car- 
rier energy is lower if they move in a FM region. Therefore, by interaction 
with the spin system in the lattice, they are able to change its magnetic 
order, creating FM micro-regions and become trapped in them. In a degen- 
erate magnetic semiconductor this is a cooperative phenomenon: a number 
of carriers are self-trapped in the same FM micro-region reducing the energy 
per carrier necessary to create it. When the carrier concentration is high 
enough, the phase-separated state turns out to be stable. It must be stressed 
that our calculation is valid only in the case of a degenerate semiconductor, 
i.e. over the range of conduction electron density in which this collective 
behaviour of the carriers appears. The micro-region contains a high number 
of conduction electrons. No individual (isolated) ferrons (magnetic polarons) 
exists over this range of concentrations. Unfortunately, there is no means of 
calculate rid accurately, simply it must be defined as the density from which 
these micro-regions contain a "high" number of conduction electrons. 

Within the Vonsovsky s-d model and using a variational procedure to 
calculate the energies of the different phases, Nagaev shows that the phase- 
separated state was the ground state for degenerate antiferromagnetic semi- 
conductors over a certain range of doping. One of the assumptions of the 
calculation is the absence of conduction electrons in the AF part of the phase- 
separated state (see the reference [||] for a review). In this article we modifed 
the calculations in |I[ to allow for the presence of conduction electrons in the 
AF part of phase-separated state and we shown that the energy minimum is 
encountered when all the conduction electrons are in the FM phase. 
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2 The Hamiltonian of the system 



Our starting point is the Hamiltonian of the generalized Vonsovsky s-d model: 

xo l a *V - \H^9- mS f +H c (1) 

where at , a^ a are the creation and destruction operators of a conduction 

electron with quasi-momentum k and spin projection a, is the operator 
of the core spins of the magnetic ions (d-spins) with number g, s aa > are the 
Pauli matrices and N is the total number of unit cells in the crystal. The 
crystal structure is assumed to be simple cubic, the d-spin magnitude being 
S. In the case of that the conduction electrons move on the magnetic ions the 
label h in the second term of equation (JJ) is equal to g. But in the case of the 

— * 

electron moving on a non-magnetic atom then h corresponds to the nearest 
magnetic neighbors of this ion and summation is carried over them. The 
term bilinear in the d-spin operators will be called the exchange Hamiltonian, 
though in reality it corresponds to the super-exchange interaction between 
magnetic ions. The term Hq describes the Coulomb energy of interaction of 
the conduction electrons with each other and with ionized impurities. The 
charge of the latter, compensating the charge of the conduction electrons, is 
assumed to be distributed uniformly over the sample (the jellium model). 

The main parameters in the Vonsovsky s-d model are W = 2zt, the car- 
riers band-width, AS the exchange energy between conduction electrons and 
magnetic ions (s-d exchange energy) and zIS 2 the exchange energy between 
magnetic ions (d-d exchange energy). / < is the exchange integral between 
first-nearest neighbors magnetic ions. The smallest parameter is the d-d ex- 
change energy, which is of the order of magnitude of Neel temperature of 
the pure compound (without doping). We differentiate two possibilities de- 
pending on the relative value of W and AS. In the case W » AS we have 
a wide-band semiconductor (e.g. EuTe). In the opposite case W « AS, 
we talk about a "double-exchange" material (e.g. lanthanium manganites). 
For the sake of definiteness the sign of A is assumed to be positive. We also 
consider only the case H = because at least up to fields H <2J (the field 
at which both sublattices collapse) the difference in energy for FM and AF 
parts always favored the conduction electron being in the FM part. 
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3 Variational procedure 



As in references QTJ, |2], |3| a variational procedure will be used here to find 
the ground state energy of the phase-separated system. Three variational 
parameters are used, the ratio x of the volumes of the antiferromagnetic and 
ferromagnetic phases, the radius R of the spheres which form the multiply 
connected phase and y = (Np — Na)/(Np + A/a), where Np (A/a) are the 
number of conduction electrons in the FM (AF) phase, y gives the ratio 
between the excess of electrons in a given phase and the total number of 
electrons and varies from —1 (all the electrons in the AF phase) to +1 (all 
the electrons in the FM phase). 

With these definitions the densities of conduction electrons in each phase 

are: 



Np 1 + y 

n F = — =n(l+x)^— 

nA = v- A =n \—)— (2) 

where n is the mean electron density in the crystal. 
The trial energy is given by the expression: 

E ps = E kin + E sur + E c + E mag (3) 

where E kin is the standard bulk kinetic energy of the conduction electrons, 
E sm is the surface electron energy, i.e. the correction to E kin that appears 
due to the quantization of the electron motion in regions of finite dimensions 
(in fact, this energy takes into account the electron level spatial quatization 
in a bounded region in the Born-Oppenheimer approximation), Ec is the 
electrostatic energy due to the inhomogeneous electronic density and E mag is 
the magnetic energy. E maug is decomposed in two parts: the energy associated 
to the s-d exchange interaction, -E s -d and the one associated to d-d exchange 
interaction, E&_&. 

We calculate now these energies for the phase-separated and the homo- 
geneous states. As discussed in the introduction, we expect two different 
topologies of the phase-separated state depending on the conduction elec- 
tron concentration. The conduction electron density at which this change 
occurs is denoted by ny. Therefore we divide our calculation for the en- 
ergy of the phase-separated state in the n d < n < n T and n T < n < n u 
cases. As we only consider here conduction electron concentrations over the 
rtd < n < n u range and magnetic fields smaller than the saturation field of 
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the material, the only possible homogeneous state we have to consider is the 
AF one. All the energies are calculated per unit volume of the unit cell and 
the numerical values are given in eV. 



3.1 Energy of the phase separated state for < n < 

Over this concentration range the phase-separated state is expected to be an 
AF matrix in which FM spheres are embedded. The total kinetic energy per 
unit volume of the unit cell is: 

E kin = + ^ (4) 

1 + x 1 + x 

where E kin (E kin ) represent the energies of the conduction electrons in the 
FM (AF) phase. 

In the ferromagnetic phase all the charge carriers are assumed to be spin- 
polarized. For W « AS this condition is met for all the charge carrier 
densities. For W » AS it is met when AS > // where fi is their Fermi 
energy. It is also necessary to take into account the different band bottoms 
encountered by the charge carriers moving in the FM and the AF part of the 
crystal ||. Therefore, in the free electron approximation, the kinetic energy 
of the conduction electrons in the FM phase is: 

= + 5^^ nF (5) 

From now on t = 2 » 2 an d M n ) = t(Qir 2 n)^. Then: 

EL = -6*n(l + + jUn)n(l + *)§ (^) ' (6) 

The kinetic energy of the conduction electrons in the antiferromagnetic part 
is: 

f -Q tn l±zl^ + p-l^) n (l±*)h±-y)l if W»AS 
E A — I x x (7) 

The surface energy, E SUT , appears due to the fact that the conduction 
electrons are confined in a finite volume. For each phase, it is proportional to 
A, being S the surface between both phases and V* the volume of each phase, 
FM, Vp or AF, Va respectively. We must point out that actually no surface 
energy independent of the conduction electron density can appear, because 
this term turns out meaningless in the case of the conduction electron density 
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equals to zero 0, where phase separation does not exist. The contribution 
to the surface energy due to magnetic interactions is also zero in the first- 
neighbour approximation considered by us. 

E F tE a 

771 sur | ^-^sur i , q\ 

E SU v = — h — — (?) 

1 + X 1 + X 

where Ef nT (E^.) represent the surface energies of the conduction electrons 
in the FM (AF) phase. 

Taking into account the results of 0, we have in the FM part: 

e f =^(iy^^L (9) 

^ sur 16 \6J V F _§ [J) 
ftp 

Using the fact that ^ = -| being R (in a units) the radius of the ferromag- 
netic spheres containing the conduction electrons, we obtain: 

1 4 

9 /7r\3 1 . . 2. , i /l + y\a . . 



For the AF part and making the obvious substitution in equation (|T0D, we 
have: 



[^(1)^2-1^1(^)1(^)1 if W»AS 

E SU P =\ ! (11) 

I ^^Yl^mH 1 ^ ^ W«AS 

The Coulomb energy Eq is calculated using the jellium model for the 
ionized impurities. The crystal is separated into the Wigner cells, i.e. into 
spheres enveloping the FM inclusions drawn so as to make the total charge 
inside the sphere vanish. This method provides a good approximation to the 
electrostatic energy for small FM volumes, i.e. in the limit x » 1. Moreover 
it supposes to admit that the crystal is isotropic and homogeneous regarding 
the spatial distribution of these spheres. So we cannot describe the fractal 
nature of the boundary that separes phases. In this approximation: 



2tt e 2 n 2 R 2 (\ - y l+y\ 2 , i, X 2 X . . 

E c = - ^ H^-*^ 2 + 3X-3X3 l + z 3 12 

5e r x A a V 2 2 J 

As was said above, the magnetic energy has two contributions. One part- 
correspond to the antiferromagnetic interaction between the ionic d-shells, 
-^d-d- We estimate it within the mean field approximation. Because our 
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calculation is constrained to T = 0, we only have to take into account the 
ground state of the system. For the antiferromagnetic phase we have: 



2 

where \J\ = For the ferromagnetic phase: 



Etd = (13) 



£ d F - d = ^ (14) 
Therefore the total energy E A _& is: 

xEt d Ef_ d \J\Sl-x 

Ed ~ d = TTx~ TTx~ = 1TTTx~ (15) 

We assume that the AF part is in an collinear AF state. A canted AF state 
would be also possible. In the fourth parameter, namely the canting 

angle, must be added to the minimization procedure. Its value depends 
only on the conduction electron density, which in the AF part of the phase- 
separated state depends on the parameters x, y. We ruled out this possibility 
in the view of the results in reference [0 where it is shown that the canted 
state is unstable against charge-carrier density fluctuations and it is never 
realized. The second part corresponds to the energy due to the interaction 
between the conduction electrons and magnetic ions, -Eg-d- Again we have: 

xE^_ d Ef_ d 
E s -d = — h — — (16) 



1 + x 1 + x 



where in the FM part: 



£ s F -d = -4^ F = -4^(i +*) l -V L (17) 



And in the AF part: 



if W » AS 
~n^^ if W«AS 



Notice that for W « AS, A does not appears in the calculation, it is only 
an additive constant. 
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3.2 Energy of the phase-separated state for < n < n 



In the case of that the density of the conduction electrons is higher than the 
percolation density nx, a phase transition is expected with a change in the 
sample topology The sample changes from being a AF matrix in which FM 
spheres are embedded to have a simply connected FM part with AF spheres 
in. To take this point into account we have to modify surface and Coulomb 
energies as follows: 



2?r e 2 n 2 R 2 x ( 1 + y l-y\ 2 fn n n i ,„ ^2 

EL = r 6 ^fli Mn)nix(1 + x)i ^¥f (19) 

A _ \ ^(f) J B2-l/i(n)nS(i±s)i(i?)J if W » AS 

I ^(f)H^(^(¥) 1 if W«AS 
The others energies remain unchanged. 



3.3 Energy of the homogeneous phase 

Over the concentration range n& < n < n u the homogeneous antiferromag- 
netic phase has an energy: 

E hom - -t%in + -^mag l 2U J 

where 



p hom 



—6tn + |2 a 2 fx(n)n 



E 



hom 



Mis 
2 



(21) 



if W >> AS and 



77 hom 



6tn 



+ ! 



fi(n)n 



V2S+1 1 5 



p hom 

mag 



L7JS 
2 



A5n 
2 



(22) 



if W« AS 



4 Results 

The stationary state of the system is determined from the condition that the 
total energy of the system equation ([3]) be a minimum with respect to the 
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three variational parameters x, R and y. The parameter R only appears in 
the surface and Coulomb energies. Minimizing the sum Q = E snT + Eq with 
respect to R for fixed x and y leads to the following expressions: 



Qmin = a>y(hfg )s 



= ° 7 (m^trJ (23) 



where: 



« = 3(f) 1 (|) § (f) f - 1-2 7=(#M a n 1 
a' = M*S^ 3 « 0.59 Y = (^^p 
We have used the following definitions: 



f(x) = 2a; + 3-3(l + 2;)3 

9(x , y) = (i±i)' (1 + ^ + x (i_i)'(l±^)» 

if nd < n < tit- 

In the case of a multiply connected AF phase, i.e. for n T < n < n u , the 
latter expressions are modified according to: 



f(x) = 2 + 3a; -3s 3(1 + a;) 3 

g(x ,y) = (l + Vj^l+.ji+^^^llf)* 

2 

In both cases we have introduced the parameter x equal to 2~s if W » AS 
and to (2S + 1)4 if << AS. 

After substituting equation (p3|) into equation (§) minimization with re- 
spect to x and y is carried out numerically. 

The percolation density is calculated by solving for the value of n at 
which both separate phases are minimal as a function of x and y and have 
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the same energy. For and n u , the requirements are that the separate state 
have a minimum as a function of x, and its energy being equal to that of the 
corresponding homogeneous state []. 

The values of the parameters that characterize wide-band semiconduc- 
tors are chosen the same as in reference || (they correspond to rare-earth 
compounds of the EuTe type): S = |, \J\S = 1CT 3 eV (that implies Tn = 5 
K, and the field at which both sublattices collapse is 98.7 kOe), AS = 1 eV, 
e r = 20, a~ 3 = 4 • 10 22 cm -3 , and the electron effective mass is equal to free 
electron mass, this implies W = 4 eV. The value obtained for = 1.10 • 10 20 
cm -3 . The value for n u , at which the conductive phase-separated state be- 
comes unstable at H = T = is 1.79 ■ 10 20 cm -3 . In this work we consider 
the following electron concentrations 0.9 tit and 1.3 n?, as an example of 
both topologies. To study the difference in energies between the energy of 
the homogeneous antiferromagnetic phase E^ om and the energy of the phase- 
separated state E ps for both electron concentrations we show in figures || 
and y a contour plot of the energy difference AE = — — hom 2 „ ps — as 
a function of x and y. Notice that AE = if E ps > -Ehom- The lines 
on the diagram are lines of constant AE. The regions with no lines corre- 
spond to AE = 0. The lines in the upper part of the figures correspond 
to more negative energy differences. It can be observed that as y increases 
AE becomes more negative. In both cases, the minimum value of AE is 
obtained for y = 1, i.e. = 0. For the values of n we show here, the 
minima are £ps ~ r f hom (x = 1.3662, y = 1) = -0.1707 at n = 0.9 n T and 
gps-ghom ( x = 0.5169, y = 1) = -0.2095 at n = 1.3 n T . 

The values of parameters that characterize "double-exchange" semicon- 
ductors are the following (typical for lanthanium manganites, for example): 
S = 2, D = -zIS 2 = 10" 2 eV (that implies T N = 58 K, and the field 
at which the two sublattices collapse at T = is 1.73 MOe), e r = 5, 
a -3 = 4 • 10 22 cm -3 , the effective electron mass equal to the free elec- 
tron mass, that implies W = 4 eV. A is regarded as high enough to meet 
W « AS, although its precise value is not needed for calculation. We 
obtain = 5.93 • 10 20 cm~ 3 and n u = 9.46 • 10 20 cm" 3 . As before we con- 
sider the following electron concentrations 0.9 tit and 1.3 n^. The energy 
minima are ( x = 1.3159^ = 1) = -0.04454 at n = 0.9n T and 

E ps -E hom ^ = 0.4681, y = 1) = -0.1162 at n = 1.3 nr. Contour plots are 
shown in figures |3] and In both cases minima are found in y = 1, i.e. all 
the conduction electrons are in the AF part of the phase-separated states. 

Analogous results are obtained for any value of the concentration of con- 
duction electrons over the range < n < n u . Therefore, we have shown that 

1 Our values for tit are slightly different from those in H, calculated imposing x = y = 1. 
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for both kinds of compounds and over the above mentioned range of concen- 
trations, the ground state corresponds to a phase-separated state with FM 
and AF parts, with absence of electrons in the antiferromagnetic part. 
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Figure 1: Countour plot of AE as a funtion of x,y for W >> AS and 
n = 0.9 nx- Lines in the upper part correspond to more negative AE. 
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Figure 2: Countour plot of AE as a funtion of x,y for W >> AS and 
n = 1.3 nx Lines in the upper part correspond to more negative AE. 
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Figure 3: Countour plot of AE as a funtion of x,y for W « AS and 
n = 0.9 nx Lines in the upper part correspond to more negative AE. 
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Figure 4: Countour plot of AE as a funtion of x,y for W << AS and 
n — 1.3 Lines in the upper part correspond to more negative AE. 
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